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Abstract 

In order to increase the accelerating gradient of Super- 
conducting Radio Frequency (SRF) cavities, Magnesium 
Diboride (MgB2) opens up hope because of its high transi- 
tion temperature and potential for low surface resistance in 
the high RF field regime. However, due to the presence of 
the small superconducting gap in the it band, the nonlinear 
response of MgB2 is potentially quite large compared to a 
single gap s-wave superconductor (SC) such as Nb. Under- 
standing the mechanisms of nonlinearity coming from the 
two-band structure of MgB2, as well as extrinsic sources, is 
an urgent requirement. A localized and strong RF magnetic 
field, created by a magnetic write head, is integrated into 
our nonlinear-Meissner-effect scanning microwave micro- 
scope [1|. MgB2 films with thickness 50 nm, fabricated 
by a hybrid physical-chemical vapor deposition technique 
on dielectric substrates, are measured at a fixed location 
and show a strongly temperature-dependent third harmonic 
response. We propose that at least two mechanisms are 
responsible for this nonlinear response, one of which in- 
volves vortex nucleation and penetration into the film. 

INTRODUCTION 

The discovery of superconductivity in MgB2 in January 
2001 [2 1 ignited enthusiasm and interest in exploring its 
material properties. Several remarkable features, for ex- 
ample a high transition temperature (T c ~ 40 K ), a high 
critical field, and a low RF surface resistance below T c , 
shows great potential in several applications such as super- 
conducting wires and magnets. The success of making high 
quality epitaxial MgB2 thin films provides another promis- 
ing application as an alternative material coating on super- 
conducting radio frequency (SRF) cavities [3|. Over the 
past decade, the improvement of accelerating gradient in 
Niobium (Nb) SRF cavities has almost reached the BCS 
limit, 57 MeV/m iHJ. In order to go further, new high T c 
materials with low RF resistance are required for interior 
coating of bulk Nb cavities. High quality MgB2 thin films 
may satisfy the demands for SRF coating materials because 
these high quality films can avoid the weak link nonlinear- 
ity between grains, and lead to the possibility of making 
high-Q cavities 0. 

However, there still exist mechanisms that produce 
non-ideal behavior at low temperatures under high RF 
magnetic fields, such as vortex nucleation and motion in 
the film (6). In addition, due to the it band and a band, 



the intrinsic nonlinear Meissner effect of MgB2 is large 
compared to other single-gap s-wave superconductors 
171 . Therefore the study of MgB2 microwave nonlinear 
response in the high frequency region (usually several 
GHz in SRF applications) can reveal the dissipative and 
nondisipative nonlinear mechanisms and allow application 
of these high quality MgB2 films as cavity coatings. 

In our experiment the localized harmonic response of su- 
perconductors is excited by a magnetic write head probe 
extracted from a commercial magnetic hard drive HI. 
Based on the gap geometry of the magnetic write head 
probe, sub micron resolution is expected. We present our 
observation of at least two measurable nonlinear mech- 
anisms involved in high quality MgB2 films below T c . 
These films were grown on sapphire substrates by hybrid 
physical-chemical vapor deposition technique (HPCVD). 
A detailed description of the growth technique has been re- 
ported before [8]. Finally, experimental nonlinearity data 
will be interpreted as a combination of intrinsic nonlinear 
response (9J and vortex nonlinearity |6|. 
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Figure 1 : Set up of phase-sensitive measurement in non- 
linear microwave microscopy. The frequency offset mode 
(FOM) of a vector network analyzer (VNA : model PNA- 
X N5242A) is used in this measurement. Inset 1. is a 
schematic picture of the probe above the sample and in- 
set 2 shows the magnetic write head probe assembly on top 
of a superconducting thin film. 



EXPERIMENTAL SETUP 

The experimental setup for amplitude and phase mea- 
surements of the superconductor harmonic response is 
shown in Fig. \T\ An excited wave (fundamental sig- 



nal) at frequency / comes from the vector network ana- 
lyzer (VNA) and is low-pass filtered to eliminate higher 
harmonics of the source signal. This fundamental tone is 
sent to the magnetic write head probe to generate a lo- 
calized RF magnetic field on the superconductor sample. 
Two insets in Fig. Q] shows close-up views of our mag- 
netic write head probe on superconducting samples. Due 
to the intense nature of this field, the superconductor re- 
sponds by generating currents at both the fundamental fre- 
quency and at harmonics of this frequency. The generated 
harmonic signal is high-pass filtered to remove the funda- 
mental signal Vf and an un-ratioed measurement of V3 / is 
performed on port 2 of the VNA. In order to get a phase- 
sensitive measurement of the 3 rd harmonic signal coming 
from the superconducting sample, a harmonic generation 
circuit is connected to provide a reference 3 rd harmonic 
signal, and the relative phase difference between the main 
circuit and reference circuit is measured. Further detail 
about this phase-sensitive measurement technique can be 
found in Ref. iflOl . In this way we measure the complex 
third harmonic voltage of y^ m P le ^ j>j or correspond- 
ing scalar power Pg mple (T). The lowest noise floor in our 
VNA is -145 dBm for the un-ratioed power measurement. 
A ratioed measurement of the complex V^f mp e (T)/V^f 
is also performed at the same time. In this paper we only 
discuss the unratioed measurements of P^ mpe (T) and 
qualitatively discuss the mechanisms of third harmonic re- 
sponse of the MgB2 film. 

THIRD ORDER NONLINEAR 
MEASUREMENT RESULTS 

The measurement of the 3 rd order harmonic power (P3 /) 
is performed near the center of the epitaxial MgB2 film of 
thickness 50 nm. The T c of this sample is 36 K as mea- 
sured by the four point resistance method. Figure [2] shows 
the temperature dependent P%f(T) curves at the excited 
frequency 5.33 GHz and excited power +14 dBm. Above 
40 K a very small signal begins to arise above the noise 
floor of the network analyzer. This P3 / is from the mag- 
netic write head probe itself. We have measured the P3 f of 
the magnetic probe on the surface of a bare sapphire sub- 
strate and in general this probe nonlinearity is negligible 
at excited powers under 14 dBm. Although excited pow- 
ers above +14 dBm excites stronger nonlinearity from the 
probe, this nonlinearity is almost temperature independent 
in the Helium cooling temperature range. Therefore probe 
nonlinearity can be treated as a constant background signal 
above the noise floor of the spectrum analyzer. The mecha- 
nism of probe nonlinearity is the hysteretic behavior of the 
yoke material [ 1 1| and has been discussed previously 1 1 1. 

From Fig. a clear P3 / (T) peak centered at 35 K shows 
up above the noise floor. This peak arises from the intrin- 
sic nonlinear Meissner Effect (NLME) at T c due to the en- 
hanced sensitivity of superconducting properties as the su- 
perfluid density decreases to near-zero levels. This peak 
at T c is also phenomenologically predicted by Ginzburg- 
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Figure 2: Temperature dependence of 3rd harmonic power 
P3/ from a 50 nm thick MgB2 measured with an excited 
frequency of 5.33 GHz. 
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Figure 3: (a) Power dependence of P 3 f on P f for the 50 nm 
thick MgB 2 film, (b) Fitted slope at selected temperatures 
for the film. The marked dip may describe a vortex-free 
region and its small slope is likely due to the probe nonlin- 
earity. 

Landau theory, and is discussed further below. 

We also note the onset of a low temperature nonlinearity 
below 27 K, which implies that another temperature de- 
pendent nonlinear mechanism is active. It may be that the 
applied RF field from the probe is strong enough to pene- 
trate into the superconductor and create deep flux penetra- 
tion or even Abrikosov vortices in a localized area. This 
new nonlinear mechanism dominates the overall measured 
nonlinearity in the low temperature region. Further quali- 
tative discussion of the low temperature nonlinearity from 
a possible Abrikosov vortex critical state will be addressed 
in detail below. 

In addition, in the temperature regime of 29K ~ 33K, 
there is a minimum P%f signal, which implies no strong 
nonlinearity mechanisms in this temperature range. The 
supercurrent of a vortex circulates around the normal core 
with an approximate size of the magnetic penetration depth 
A(T). Therefore once X(T) is bigger than the film thick- 
ness (50 nm in this case) above a certain temperature, vor- 
tex penetration due to parallel magnetic field will be sup- 
pressed fl2l . This regime would be very suitable to fab- 
ricate a low nonlinearity superconducting response in a 
multi-layer superconductor / insulator structure ifPJl . 

Measurements of the dependence of P3/ on Pf are 



shown in Fig. [3] (a) for the 50 nm thick MgB2 film at 
some selected temperatures. In the normal state of MgB2, 
the measured nonlinearity comes from the probe itself and 
shows a slope steeper than 3 at high excited power above 
+ 15 dBm. In the intrinsic nonlinear Meissner regime, the 
slope is 2.84, very close to 3 as predicted for the intrinsic 
NLME 1 14 1 . Based in part on this evidence, we believe that 
in the high temperature region close to T c , the P 3 j comes 
from the intrinsic NLME. In the low temperature regime, 
the slopes of P3/ vs. Pf are around 1.5. This value is sim- 
ilar to that predicted by many phenomenological models 
(between 1 ~ 2) in an Abrikosov vortex critical state [15] 
|[T6l . It should be noted that the low temperature nonlinear- 
ity can be easily excited at low power. Figure[3](b) show the 
P3/-P/ slope evolution from a flux/vortex dominated non- 
linear regime at low temperature to a NLME regime around 
T c for this MgB 2 film. 

INTRINSIC NLME OF MGB 2 

The intrinsic nonlinearity comes from the backflow of 
excited quasiparticles in a current-carrying superconduc- 
tor, which results in an effective decrease of the superfluid 
density. Therefore, a two band quasiparticle backflow cal- 
culation should be applied to the MgB 2 intrinsic nonlinear- 
ity. Based on the work of Dahm and Scalapino[9|, the tem- 
perature and induced current density dependent superfluid 
density n s (T, J) can be written as 



n s(T, J) J 2 

it n\ = 1 - W > ' Jnl = 

n s (T, 0) J NL 



Jr. 



b^(T) + b a (T)- 



(1) 

where b CT and are the temperature dependent nonlinear 
coefficients for the a band and it band, respectively, and 
their values are defined in reference 0. J C CT = 4.87 x 
l0 8 A/cm 2 and J C)7r = 3.32 x 10 8 A/cm 2 are the pair- 
breaking current densities for the two bands. For a 50 nm 
thick MgB 2 thin film, the generated third harmonic power 
P 3 f (T) is estimated by substituting Jjy j, into the following 
equation |[T4ll 



^^(T)T 2 



(2) 



32Z d*J% L (T) 

where oj is the angular frequency of the incident wave, d is 
the film thickness, A(T) is the temperature dependent mag- 
netic penetration depth, Z is the characteristic impedance 
of the transmission line in the microscope, and T is ge- 
ometry factor which is estimated to be 10 5 ^4 3 /m 2 for the 
magnetic write head field distribution under a 100 mW ex- 
cited power. The solid red line in Fig. [4] shows the P%f{T) 
simulated results of Eqs. (Q} and (O for the 50 nm thick 
film at a 5.33 GHz excited frequency. This intrinsic NLME 
response has measurable values above the noise floor only 
in the high temperature region near T c . The experimental 
data of the MgB 2 film under a +18 dBm, 5.33 GHz mi- 
crowave excitation is shown in the blue dots. Therefore, at 



-100 

^-110 
a 

£9-120 

^-130 
-140 
-150 




5 10 15 20 25 30 35 40 
Temperature (K) 

Figure 4: Solid red line is the simulated result of the in- 
trinsic NLME of MgB 2 with thickness 50 nm under the as- 
sumption that the magnetic probe provides a field described 
by a geometry factor T = 8.3 * 10 5 A 3 /m 2 . Other pa- 
rameters used in this calculation include A(0_ft')=100 nm, 
A ctlto// =800nm, J cutof f = 4.2*10 n A/m 2 andT c =34.6 
K with a standard deviation of Gaussian spread of 5T C =1.3 
K. The blue dots are the data from the experiment. 



lower temperatures the nonlinear mechanism must be of a 
different nature. 

NONLINEARITY IN ABRIKOSOV 
VORTEX CRITICAL STATE 

Nonlinearity From Moving Vortices 

Vortex nucleation and penetration into the film induces a 
dynamic instability and generates harmonic response. The 
equation of motion of a vortex in a semi-infinite supercon- 
ductor driven by a harmonic magnetic field is given by 1 6 1 



MoA(T) 



27174, A 3 (T) 1V A(T) 



where x is the coordinate of the vortex position with respect 
to the surface (x=Q), n is the Bardeen-Stephen vortex vis- 
cosity, $0 is me fl ux quantum, uj is the angular frequency 
of the incident wave, no is the permeability of vacuum, Bq 
is the magnitude of RF magnetic field on the SC surface 
and K\ (x) is the modified Bessel function. The first term 
on the right hand side is the Lorentz force per unit length 
on the vortex due to the screening currents created by the 
driving field. The second term on the right is the force per 
unit length exerted by the image vortex that arises from the 
SC/vacuum surface. This equation assumes a bulk super- 
conductor. 

The solution for the trajectory of this single vortex is 
shown in Fig. [5] as a function of time in the lower solid 
blue curve. The applied RF field B(t) is also included in 
the figure to illustrate the relation of the vortex position and 
the applied field with time. The time for the first vortex en- 
try can be determined as 0, 



tn = 



arcsin(|^) 



(4) 




Figure 5: A calculation of vortex position x(t) with the 
applied RF field B(t) over several RF cycles. Here to in- 
dicates the time of vortex entry, t m is the time when the 
vortex reaches its maximum distance into the film, t c is the 
antivortex nucleation time, and t a is the vortex / anti-vortex 
collision time. 

where B v is the penetration field of a vortex (assuming B v 
< Bo). A vortex will start to nucleate and enter into the 
film when B(t) exceeds the Bean-Livingston barrier ifTTl . 
This vortex also creates a supercurrent circulating around 
the core and distorts the Meissner screening current near 
the surface. During the reverse part of the RF cycle, the 
Meissner screening current is enhanced so that at time t c an 
anti-vortex will penetrate into the superconductor as shown 
by the dashed red line. This second vortex will annihilate 
with the first vortex at time t a . This procedure of vortex- 
antivortex entry and annihilation continues and will gener- 
ate a third harmonic signal. 

For further quantitative modeling, the following two 
cases should been taken into consideration: 

(1) Finite thickness of the film. 

The vortex equation of motion given above is only suit- 
able for bulk materials and assumes that a uniform mag- 
netic field is applied parallel to the SC surface. In the finite 
thickness case, an infinite number of image vortices are re- 
quired to satisfy the boundary conditions. However we can 
make an approximation that only two image vortices are 
required. Therefore the equation of motion of the vortex is 
modified to, 

/^oA A 

_^ i( 2^_*) }] (5) 

where d is the thickness of the film. This modification for 
a second image force will help improve the quantitative 
modeling. 

(2) surface roughness of the film 

In the Bean-Livingston model the superconducitng sur- 
face is assumed to be a perfect plane IfTTl . When the 
surface has roughness with characteristic length > £ 
(coherence length), a geometry effect should be taken 



into consideration [18]. Generally, for a sharper corner, 
the Meissner screening current density will be enhanced 
and the penetration field of the first vortex entry {B v ) will 
decrease. For example, at a corner with a 90° angle, an 
enhancement of the screening current is roughly estimated 
to be a factor of 4 IfTSl . This means a vortex will penetrate 
at sharp points or cusps easily and reduce the vortex 
nucleation time during the RF cycle. Therefore, nonlinear 
harmonic response will be increased compared to the case 
of a perfect plane. Hence for a given excitation level, the 
harmonic response will depend on the surface topography, 
and an image showing this contrast can be built up by 
raster scanning the magnetic probe. 

Nonlinearity From Switching Between the Meiss- 
ner State and the Vortex Critical State 

In addition to vortex and antivortex nucleation and mo- 
tion, another possibility to generate a P^f in the Abrikosov 
vortex critical state is the switching between this state and 
the nonlinear Meissner state. While the peak value of the 
applied RF magnetic field is higher than the surface pene- 
tration field of the superconductor, the material will switch 
into the critical vortex state from the Meissner state. This 
process of switching between states implies another source 
of nonlinear harmonic response. 

Fig. [6] (a) shows a schematic illustration of our exper- 
iment in which the RF magnetic field from the magnetic 
write head probe interacts with the superconductor under- 
neath the probe. One can model the flux distribution with 
an equivalent magnetic circuit as shown in Fig|6](b). The 
inductively coupled driving line provides a magnetomotive 
force (V m ) to the yoke with a reluctance R y . A magnetic 
flux $ is channeled down along the yoke to the gap. There 
the flux can divide into two branches: one directly goes 
through the gap with a reluctance R g and the other shunts 
into the superconductor with a reluctance R sc . The reluc- 
tance of superconductor R sc is a time-variable reluctance. 
It is a combination of the reluctance R s from the nonlinear 
Meissner state and the reluctance R v from the vortex crit- 
ical state. While the applied field B(t) is smaller than the 
penetration field B v , the reluctance will remain at the value 
of R s . Once B(t) > B v , an additional reluctance channel 
R v is created. Whether the vortex enters as a semi-loop (as 
assumed above), or as a vortex- antivortex pair, remains to 
be evaluated. 

Because a magnetic circuit is analogous to an electric 
circuit, we can compute R v and R s with node-voltage anal- 
ysis. Assume that the flux going through R g and R sc is <p g 
and c/) sc , respectively. In the nonlinear Meissner state, we 
obtain 

(j) g = B Q A gap = Bo{w * g) ; (f> sc = B (w * A) (6) 

where A gap is the cross-sectional area of the gap, w and 
g are the width and the thickness of the gap, respectively 
(see Fig.|6](a)), and A is temperature dependent penetration 
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Figure 6: (a) Schematic illustration of the magnetic flux 
coming from the yoke to the superconductor where w,g 
and l g represent the width, the thickness, and the length of 
the gap, respectively. The length l g is on the order of 200 
nm for our write head probe, (b) An equivalent magnetic 
circuit for the magnetic flux transport from the yoke to the 
superconductor. Note R SC =R S +R V is a time variable reluc- 
tance.. 



depth. Applying the node-voltage law, we have (f> g R g = 
(j> sc Rs- Finally, the reluctance of the nonlinear Meissner 
state is given by 



R s = -rR g 



Rg = 



In 



H Ag 



(7) 



where l g is the length of the gap (see Fig. [6] (a)) and fi 
is the permeability of vacuum. On the other hand, in the 
vortex critical state, we have the magnetomotive force as 



t>oR v — I 

J na 



H Q dy 



Bq1 9 
Ho 



BqX (8) 



where (f>o is the flux quantum and the latter equation as- 
sumes that the vortex carries one unit of magnetic flux, 
which may not always be the case. Hence the reluctance 
of the vortex critical state is R v = — H-r. In addition to 
the nonlinearity of the penetration depth with RF field, we 
believe the transient between R v and R s will also induce a 
third harmonic response. 



CONCLUSIONS 

A strongly temperature-dependent third harmonic re- 
sponse is found in high quality MgB2 films. In addition 
to the intrinsic nonlinearity, the nonlinearity coming from 
the Abrikosov critical state may also be involved. From 
the dependence of P3/ onPj, the nonlinearity mechanism 



changes from a intrinsic nonlinear Meissner effect to a pos- 
sible vortex critical state dominated nonlinearity upon cool- 
ing the high quality epitaxial MgB2 film. The mechanics 
of nonlinearity in the Abrikosov vortex critical state can 
be qualitatively interpreted by two models - first: annihila- 
tion of moving vortex & antivortex pairs and second: state 
switching between a Meissner state and a vortex critical 
state. 
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